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ABSTRACT 

Recent observations reveal that magnetic turbulence in the nearly colisionless solar wind plasma 
extends to scales smaller than the plasma microscales, such as ion gyroradius and ion inertial length. 
Measured breaks in the spectra of magnetic and density fluctuations at high frequencies are thought 
to be related to the transition from large-scale hydromagnetic to small-scale kinetic turbulence. The 
scales of such transitions and the responsible physical mechanisms are not well understood however. 
In the present work we emphasize the crucial role of the plasma parameters in the transition to kinetic 
turbulence, such as the ion and electron plasma beta, the electron to ion temperature ratio, the degree 
of obliquity of turbulent fluctuations. We then propose an explanation for the spectral breaks reported 
in recent observations. 

Subject headings: magnetic fields — magnetohydrodynamics — turbulence 


1. INTRODUCTION 

In situ measurements of the magnetic, electric, and 
density fluctuations in the solar wind provide valuable 
information on nonlinear dynamics of a nearly colli¬ 
sionless astrophysical plasma. At large hydrodynamic 
scales (corresponding to low frequencies in the measure¬ 
ments), such fluctuations are thought to be consistent 
with magnetohydrodynamic turbulence viewed as inter¬ 
acting oblique Alfven modes propagating along the back¬ 
ground magnetic field (IIroshnikovlll963tlKraichnanlll965l: 


Goldreich fc Sridhail 119951 IGaltier et al.ll200Ct iBoldvrevI 

200ti[l . At higher frequencies the spectrum of such turbu¬ 


lence exhibits a break, which corresponds to the spatial 
scale broadly consistent with the plasma microscales such 
as the ion gyroradius or the ion inertial length. 

It has been proposed that the spectral break 
in the solar wind and other astrophysical plas¬ 
mas can mark a transition from the non-dispersive 
Alfven modes to the dispersive kinetic-Alfven 


modes (iBale et al.l 120051 ILeamon et al.l I1998L 

1999; 

HollweellI999l iHowes et al.l 120061 Chandran et al.l 

2009; 

Shaikh & ZankI 20091: Ghandran et al.ll20I0l: Chen et al. 

20IC 

; IHowes & Ouataert 120101: Petrosvan et al. 

20IC 

; IHowes et al.l 120111: Chandran et al. 

2011; 

TenBaree & Bowed 120121 IBoldvrev & Pered 

2012; 

Sahraoui et al. 120121 IMithaiwala et al.l 

2012; 

Boldvrev & PerezI 120131: iPodestal 120131: IChen et al. 

201.31: iHaverkorn & Soansfled |201.3(1. A possibility of 


transition to whistler turbulen c e has also been con- 
sidered (iBeinroth &: Neubaued 119811: IGoroniti et al.l 


1982t IGoldstein et al. 


IStawicki et al.l 


Galtier fc Bhattacharied I2003L l2005t Gary et al.l 


2001 


2008 


Saito et al.l 120081: Gary &: Smith! I2009t iShaikhI 2009 . 
20101: ^rv et al.[l2010D. howeve r, recent studies (e.g., 
Podestal l2013t iGhen et al.1 1201, 3ll suggest that whistler 
turbulence, if present at subpropton scales, contributes 
only a small fraction of fluctuations energy. A conclusion 
is then drawn that the spectral break occurs at the 
proton gyroscale. 

A recent work by iChen et all (1201411 tested this pre¬ 


diction by analyzing the solar wind intervals having very 
large and very small plasma beta, the ratio of the ki¬ 
netic energy of the plasma particles to the magnetic en¬ 
ergy. In particular, the intervals were selected with ion 
and electron plasma beta satisfying ^ /3e ^ 1 and 
1 ^ /3e ^ A- In the first case, the theory of oblique 
Alfvenic turbulence predicts the break at the ion gy¬ 
roscale, in the s econd one at the ion acoustic scale. The 
observations of (jGhen et al.ll2ni4 l agree with the theory 
in the first case, and disagree in the second, where the 
break is observed at the ion inertial length instead. This 
puzzling result may question the applicability of the the¬ 
ory of turbulent cascade to the solar wind plasma. 

We address this contradiction by inspecting the theory 
of Alfven turbulence in the limiting cases of large and 
small plasma beta. We consider various mechanisms 
that may be responsible for the spectral break, includ¬ 
ing the possibility that the major assumption of the 
standard theory, the obliquity of propagation, can break 
down in the case 1 » /3e 2> /3i. The reason for the latter 
possibility is upscatter of the Alfvenic fluctuations due 
to their interactions with the ion-acoustic modes and the 
fast modes that are weakly damped in a non-isothermal 
low-beta plasma. The Afvenic turbulence then develops 
a non-oblique component fc|| > k±, which is dissipated 
due to the ion cyclotron resonance at the ion inertial 
scale thus explaining the observed spectral break in this 
case. 


2. KINETIC DERIVATION 

In what follows we will use the notation: uipa = 
1^/4717100< 7 ^/mo, is the plasma frequency, and VTa = 
nia is the thermal velocity associated with the par¬ 
ticles of kind a. For a plasma consisting of electrons and 
ions, the so-called ion-acoustic velocity can be dehned. 
Vs = It is also convenient to introduce the 

Alfven speed va = Bq/ y/Amtomi, and the plasma beta, 
which is the ratio of the thermal energy of the particles 
to the magnetic energy of the plasma, and which can be 
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different for the ions and the electrons, Pi = 

Pe = 2v1lv\. The particle gyroradii are denoted by 
Pa = VTa/^a, where fla are corresponding gyrofrequen- 
cies; the ion inertial length is denoted by di = VA/^i- 
In order to find the dispersion relations for the waves 
that can propagate in a plasma, one needs to solve 
the equation det{Dim) = 0, where Dim = k^^im — 
hkm — ^e/m(w,k), and the general expression for di- 
electric tensor eim can be found in (e.g.. iStixl 119921 : 
[Aleksandrov et al1ll984 iBoldvrev et al.ll2013f) . Analytic 
solutions of this equation can be found for the important 
limiting cases that we consider below in more detail. 


2.1. Spectral breaks at Pi ^ 1 

In this case the dispersion corrections associated with 
the ion gyroradius, ~ kpi dominate the corrections asso¬ 
ciated with the ion inertial length, ~ kdi, and the lat¬ 
ter can therefore be neglected. In terms of frequency of 
Alfven waves this means that w <C fli and w <C k^VTi- 
Keeping in mind that we will be interested in oblique 
propagation, we will also assume k^Pi <C 1, however, we 
do not make assumptions about smallness of k±pi. 

The condition a; <C Hi allows one to simplify the dielec¬ 
tric tensor. The approximate expressions for the compo¬ 
nents of Cij including both electron and ion contributions, 
to the first order in small k^p^ have the form: 


k^ 7,2 


1 - Ao{z±) + &Z\\ 


A„(z l) 


n—1 


^yy ~ ^xx 


4iUjp^z± A'^(z±) 




n—1 






^XV - ^ 


^XZ — ^zx 


‘^^pikz An(z±) 


uj^,k± 




wflikz 


[1 -I- Aq(2;x)] 


I TT 

+ \lz: 


^pik± 


2 ^ikz\kz\vTi 
UJ. 

z'^Ti 




>e , '^pi n \ An{z±) 

+ IF A 


klvTe 


^ n—1 


71 UJ 


pi 


2 \kz\vTi klv"^ 


Ao{z±). 


Ti 


( 1 ) 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


( 6 ) 


In these expressions the coordinate frame is chosen such 
that the z-direction is along the uniform magnetic field, 
and the wavevector has the coordinates {k±,0,kz)- We 
use the notation z± = k^pf, Z|| = k^p^, A„(zj_) = 
In{z_\_) exp(— 2 j_), where /„ is a modified Bessel function 
of order n, and A'„{z±) denotes the derivative with re¬ 
spect to z±. The last terms in expressions @, ®, and 
® describe the collisionless dissipation by ions. As we 
will see momentarily, the imaginary term in Cyy becomes 
significant close to the ion gyroscale. 


The Bessel functions in the dispersion equation 
det(II;m) = 0 can be simplihed for large and small argu¬ 
ments. We consider these cases separately. In the case 
z_L <C 1, one can check that for oj k^VA and Pi 1, 
the dissipation term provides the dominant contribution 
in both Cyy and Dyy, while it can be neglected in ®. In 
the same limit, the eyz,zy ^xz,zz components of the 
dielectric tensor can be neglected to the leading order 
in XjPi. 

The dispersion equation then takes the form: 

^2 ^ k\c^ ^ 

^zz ^xx ^xx^yy 


where the last term describes the dissipation correction, 
and one can substitute uj = k^VA in this term. As a 
result we obtain: 


z 1 2 2 
UJ =k^VA 


1-k 


T, 


T.+Te 


klp^ - dktpt- 


2 „2 


■ 9 1,2 2 


( 8 ) 


which is the dispersion relation for Alfven waves in a 
collisionless high-beta plasma. This relation holds for 
the shear-Alfven modes, where the inertial corrections (~ 
kdi) are negligibly small. The pseudo-Alfven modes are 
strongly damped in the considered limit, and they are not 
discussed here. We note that the dissipation rate, given 
by the expression y/a; = — {9/16)i/PiJn k'^p'f, becomes 
significant at the ion gyroscale. 

For z_L ^ 1, the non-dissipative components of the 
tensor Dim agree to leading order with the components 
of Dim, for the subproton kinetic-Alfven waves obtained 
in (iBoldvrey et al.l 120131 Eqs. (32-37^. The dominant 
ion dissipation term comes in this limit through the last 
term of ([2]). We thus get the dispersion relation for the 
high-beta kinetic-Alfven waves: 


2 klk']y\ 


1 - 


^±Pi 


( 9 ) 


The dissipation rate estimated from this expression, 
y/w = —l/(4fc^p^), also becomes significant at the ion 
gyroscale. The resulting dispersion curve is schematically 
shown in Fig. O- The dispersion and dissipation terms 
in ® in]) suggest that the spectral break occurs close to 
the ion gyroscale pi independently of the obliquity of 
propagation^ This is consistent with the observational 
measurements (jChen et al.ll20l4) . 


2.2. Spectral breaks at Pi 1 

In this case the ion inertial length is much larger than 
the ion gyroradius, and we can neglect the dispersion 
corrections associated with the ion gyroradius (the finite 
Larmor radius effects). We thus assume kj^pi <C 1, but 
for the frequencies we now have kzVTi ^ uj k^vre- 
Other than that our consideration at this stage is gen¬ 
eral in that we do not require Pe, kdi, or fcz/fcj. to be 

^ The significance of the collisionless ion dissipation in the 
vicinity of the ion gyroradius has been observ ed in numerical so¬ 
lution of the Vlasov-Ma xwell equations (e^. . I Howes et ahl 120061 : 
ISchekochihin et al.l 120091 ). Our formulae ([Sl I9)l provide analytic 
explanation for those observations. 
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Figure 1. Sketch of dispersion relations for oblique waves in a 
collisionless high-/3i plasma. The dashed lines indicate the strongly 
damped modes. 


small. Our consideration here is different from previ- 
ous studies of low-be t a Alfven dispe rsion relations (e.g., 
ILvsak fc Lotkol [T996I : iHollwed 1199911 . where the hydro- 
dynamic approximation w <C Hi was made, and the ion 
inertial corrections proportional to kdi were therefore ne¬ 
glected. The dielectric tensor now has the form: 


^XX 


^yy 


UJ 


2 

pi 


ijJ‘ 


-nr 


UJ 


2 

pi 


-nv 


UJUJpi 

^XZ — ^ZX ~ 0; 


— Czy ^ 

. A 
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^pe. 


1-h 

iklv^d 

' Pz^Te 

a;2 

J 


( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 

(14) 

(15) 


where in ezz we kept the correction provided by the ion 
temperature, which will be needed for the dispersion of 
the ion-acoustic modes in section 4, but can be neglected 
for now. As a result we obtain the general dispersion 
relation for a low-,0i plasma: 


(..^ - klvl) [(^^ - k\l) - klvl) - u:W^eA] 
-u^kivl - klv\) +w^kivlk^v\kldl = 0(16) 


As noted in ( 

Kuvshinovlll994 lOhsaki & Mahaianll2004l 

IHirose et al.l 

200411. equation (|16p can also be derived 


from the Hall MHD with cold ions, Ti —>■ 0, however, 
in order to keep the finite ion-temperature effects, one 
needs to use the kinetic derivation. The terms contain¬ 
ing di describe the dispersion related to the ion inertial 
length. 

The following two limiting forms are often discussed in 
the literature. The first is the hydrodynamic limit uj <C 
fli- In this case the terms containing di can be neglected, 
and we recover the well known dispersion relations for the 
Alfven, fast and slow modes. This limit is however not 
of importance for us since it neglects the dispersion. 

The second limit is the low electron beta case /3e "C 1. 
In this limit we can neglect the terms containing the ion 
acoustic speed, that is, the second line in (ITCl) . As a 
result, the low-frequency ion acoustic mode = k^v^ 


decouples^ and from the remaining equation we can de¬ 
rive the dispersive Alfven and fast modes. 

Here we consider the case of oblique propagation 
k1/k\ <C 1. We do not make any assumptions about 
Pe though. We start with the case when the propagation 
is oblique enough so that kl/k^ < jde- (If the oppo¬ 
site inequality holds we are back to the previous limiting 
case of small /3e, we will discuss this case later.) In this 
limit the fast mode has high frequency, and it decouples 
from the other two modes. The equation for the fast 
mode is obtained if we approximate — k^v\ ~ and 
oj"^ — klvl Ri in (flBl) . which gives: 

uj'^ - u? [kiv\{l + fc^d^) -h k\_vl] -t klk\^vlv\dl = 1(17) 

For kzdi <C 1 -I- {vs/vaY the solution is the fast or mag- 
netosonic mode: 


u? = k]_{v\ + v'l), (18) 

while in the opposite case kzdi ^ 1 + {vs/vaY if trans¬ 
forms into the whistler mode: 


= klklv\dl (19) 

The remaining dispersion equation for the two low- 
frequency modes is obtained if we approximate — 
k‘^v\ Ri —k‘^v\ in (I16p . which now turns into: 

[l + kldj+vl/v\]- 

- up- [2klv‘^^ + klv\ + klvlk'idj] + k^v'^^v\ = 0 . ( 20 ) 


In the long wave limit, when the k^d^ containing term is 
negligible, the solutions are the shear-Alfven mode and 
the ion-acoustic mode. At large the shear-Alfven 

modes turns into the kinetic-Alfven mode, 


w 


2 


klk\v\p 


2 

s 



while the ion-acoustic mode becomes: 


( 21 ) 


= VlI cos^ d, 


( 22 ) 


where 9 is the angle between the wavevector and the 
background magnetic field. For each of the modes, the 
transition scale between the their long wave and short 
wave asymptotic solutions can be defined as the scale 
where the two solutions formally match. We thus ob¬ 
tain the transition scale for both the Alfven and the ion- 
acoustic modes: 


fcid?=i+4- 

vi 


(23) 


A similar result can be obtained in the framework of 
Hall MHD, e.g., (iSchekochihin et 3111200911 . The general 
behavior of these dispersion curves is schematically pre¬ 
sented in Fig. ©. 

If the propagation is not very oblique so that k\/k\ > 
/3e, then, as mentioned above, the dispersion relations 
can be found from the second limiting case. They have 


^ The ion-acoustic mode is Landau damped unless Te > 6Ti, 
which limits /3e from below. 
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Figure 2. Sketch of dispersion relations for oblique waves in a 
collisionless low-/3i plasma. The case of relatively high obliquity 
k^lk^ <C /3e < 1 is shown. The boundaries between different 
asymptotic regimes are marked by dots on the curves. 


Oblique waves in a low Pi plasma 



Figure 3. Sketch of dispersion relations for oblique waves in a 
collisionless low-/3i plasma. The case of relatively low obliquity 
< kl/k\ < 1 is shown. The boundaries between different 
asymptotic regimes are marked by dots on the curves. 

especially simple forms in the case <C 1, where we 

get for the fast and Alfven modes, respectively 

u:^ = kW{l + kldl), (24) 

uj^ = klv\l(l + kldi). (25) 

These dispersion relations are shown in Fig. ([3]). In both 
cases the breaks in the dispersion curves appear at kz ~ 
1/di- 

The obtained results suggest that for oblique propa¬ 
gation kz fcj_, and small electron beta 1 > /3e ^ ft, 
the spectral break should occur at the ion acoustic scale 
Ps^ which is smaller than the i on inertial scale ft . This 
contradicts the observations of (jChen et al.ll20l4l , where 
the break is observed at the ion inertial scale for such 
plasma parameters. In the following sections we analyze 
the physical processes that may explain this contradic¬ 
tion. 

3. TRANSITION FROM ALFVEN TO KINETIC-ALFVEN 
TURBULENCE IN THE STRONGLY NONLINEAR 
REGIME 


So far our consideration was based on the linear anal¬ 
ysis. As a justification for that we note that in the case 
of oblique propagation, the turbulence in the solar wind 
is expected to be critically balanced. This means that 
the linear frequencies of interacting modes are compara¬ 
ble with the rates of nonlinear interaction at all scales. 
Such assumption is based on phenomenological con¬ 
sider a tions (e.g . iKraichnanI 1196^ iGoldreich fc Sridharl 
Il995t iBoldvre-^ 120061 : iPerez et al.l 1201211 . and it is con¬ 
sistent with numerical studies of both Alfven and 
kinetic-Alfven turbulenc e (e.g., iCho fc Vishniad 120001 : 
iTenBarge fc Ho\^^l2012f ). Qualitative properties of lin¬ 
ear modes can therefore be present in such turbulence 
even in the strongly nonlinear regime (iB ale et al l 120051 : 
iHe et all [Ml: iHowes et al.l IMl iKlein et al l 1201211 . It 
may therefore be reasonable to relate the spectral breaks 
in the Alfvenic turbulence to the change in the dispersion 
properties of the associated linear modes. 

In the limit 1 ^ /3e ^ ft the nonlinear dynamics of 
oblique kinetic-Alfven waves is captured by the system 
of three equations for the electron density n, and the 
magnetic and velocity potentials, b i = z x V/;, an d 
vj_ = z X V<j) (|Camargo et al.l 119961 : iTerrv et 1111120011) . 
Let us normalize the independent variables as follows, 
x' = x/L, t' = tvAjL-, and the fields as ft = 4’I{Lva)^ 
n' = Sn/riQ, ip' = ^/(SqA), where L is the outer scale. 
In the dimensionless variables the system takes the form 
(we omit the ' signs): 

ftV' = V||</-^fV||n, (26) 

ftn-Fzx V/>-Vn= f-ViiViV', (27) 

X V(/-= V||Vift, (28) 

where we have denoted: Vy = dz + z x V_l^/ • V_l. 

It is seen (by rescaling n —>■ n/ftft/2, for example) 
that the only relevant microscale in the dynamics is 
ft-ft/3e/2 = Ps, the ion-acoustic scale. At larger scales, 
k±Ps “C 1, the equations for ft and ft turn into the re¬ 
duced MHD equations that describe shear-Alfven waves. 
The density decouples, and it is advected as a passive 
scalar. In the limit of small scales k±ps ^ 1, the ve¬ 
locity field decouples while density becomes dynamically 
important, and the result ing two-field system des cribes 
the kinetic-Alfven modes (iBoldvrev fc Per^l2012ll . 

To analyze the turbulence of Alfven waves we solve the 
above equations I|261I28I) in a periodic 3D rectangular do¬ 
main using a pseudo-spectral c ode. The numerical pro - 
cedure is analogous to that in (|Boldvrev fc Per^l2012f) . 
The background magnetic field is along the z axis. We 
supplement each of the equations with the random driv¬ 
ing forces fip^fn, and that are applied at small k in 
the Fourier space, satisfying 2njL± < kx,y < 47r/Lj_ and 
^'xjLz < kz < AxTjLz. The Fourier components of the 
forces are drawn from a Gaussian distribution with zero 
mean and the variance chosen as to maintain same av¬ 
erage amplitudes of the velocity, density and magnetic 
fields at the driving scale (the situation similar to ob¬ 
servations). We choose vaL± ^ UmisTj, so the turbu¬ 
lence is critically balanced at the outer scale. The aver¬ 
age force refreshing time r = where tA = LzjvA- 

The number of grid points is the same in all directions, 

Af3 = 10243 . 

The system keeps evolving until it reaches a steady 
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Figure 4. The Fourier spectra of velocity, magnetic, and density 
fluctuations in Alfvenic turbulence driven at large scales, in a low- 
beta plasma. The inertial scale corresponds to kj_ = 10, the ion 
acoustic scale is k± ^ 30. The transition from shear-Alfven to 
kinetic-Alfven turbulence is observed at the ion-acoustic scale. 

state. We select the data from over where is 

Alfven crossing time. The interval of data sampling is 
^ 0.2tA- The resulting spectra of strong oblique turbu¬ 
lence are shown in Fig. |31 It shows the spectral steep¬ 
ening indicating the transition from Alfven to kinetic- 
Alfven turbulence at the ion-acoustic scale k± Ri 30. 
The density enhancement associated with increased com¬ 
pressibility of kinetic-Alfven modes, and the “knee” in 
the density spectrum are also observed at the same scale. 
The observed break in the spectrum of strongly non¬ 
linear turbulence is less consistent with the ion-inertial 
scale, kpcrp = 10. 

4. ANOMALOUS RESISTIVITY 

In this section we discuss whether anomalous dissipa¬ 
tion may be responsible for the spectral break in Alfvenic 
turbulence. In the case of low beta and highly non- 
isothermal plasma, /3i /3e < 1, one may encounter a 

situation when the break in the spectrum is related to the 
nonlinear effect of anomalous resistivity. The physics of 
the p henomenon is the following (see, e.g.. lPaDadoDOulosl 
[197^ . As the scale of the Alfven fluctuations decreases, 
the field-parallel electron current increases, since Sb\ de¬ 
creases slower than A in a turbulent spectrum. When 
the electron velocity associated with the current exceeds 
the phase velocity of the ion acoustic waves, these waves 
become unstable due to Cherenkov radiation. The inter¬ 
action of the current with the waves then leads to the 
current damping, leading to the enhanced or anomalous 
resistivity. 

To estimate when this happens, one needs to retain 
the imaginary parts in the dielectric tensor for the ion 
acoustic waves. In a low-beta plasma, for the ion-acoustic 
waves we have uj <g; k^VA and w k±_v a- In this case, 
one can check that the and components of the 
dielectric tensor provide dominant contributions 0 . The 
imaginary parts are small in Cxx] they should be retained 

^ Another way to understand this fact is to note that the ion- 
acoustic mode is a potential mode, therefore, its dispersion rela¬ 
tion can be derived from the Poisson equation ki^ijkj = 0, which 
transforms in the adopted reference frame, k = (kj_,0^ kz), to the 
equation k'^Exx + k'^Ezz = 0. 


only in the Czz component: 





(w - kzVopl^ 

\k.\^vk 


\k, 





2k^vl^ 



where vq is the velocity of the electron current. For sim¬ 
plicity it is assumed to be constant 0 The dispersion re¬ 
lation for the ion-acoustic wave with the ion-temperature 
corrections in dTSl) has the form: 


2 

UJ^ = * ® 


l + klp2 _ 


1 + h k±Ps) 


(30) 


where the second term in the square brackets is a small 
correction. 

It is easy to see from (|29)) that for sufficiently large vo, 
the imaginary part of e^z becomes negative, leading to 
the instability. We consider separately the cases when 
the unstable modes have the wave lengths fcj_ps <C 1 and 
k±Ps ^ 1 Icf. lPaT)adoT)OulosllI977ll . In the first case the 
instability condition reads 


Vo > Vs 






(31) 


This condition holds in the case of not very strong tem¬ 
perature difference, T^/Ti < 12 for the hydrogen plasma 
(see below). In the second case the instability condition 
is 


I'D 


VTi 

k±Pi 






The right hand side achieves minimum when 


1 

(fc_L/Oj)2 



(33) 


which together with ko_Ps > 1 implies for the hydrogen 
plasma Te/Ti > 12. The instability condition is then 


Vo > VTz 




In 



-3, 


(34) 


which for T^jTi > 12 gives vo > S.bvn- 
The electron velocity vo can be estimated from the 
field-parallel electron current Jy = noevo- The current is 
related to the magnetic field fluctuations: zk x = 
(47r/c)Jfe. The rms current can therefore be found 
from the measured spectrum of magnetic fluctuations 
= (c/Any J k'^E{k)dk. Such a c urrent can be est i- 
mated from the observational sets of (jChen et al.ll20l4) . 
see Fig. [5] The estimates lead to the electron velocity vq 
about 4 to 15 times smaller then the ion-acoustic veloc¬ 
ity for the measured plasma parameters. The current is 
however known to be highly intermittent in both MHD 
and k inetic turbulence (|Zhdankin et al.l[2013t iWan et HI 
[2013), meaning that a significant fraction of the energy 
is dissipated in the regions where the current is several 


Which may be a good approximation if the excited ion-acoustic 
modes have wavenumbers and frequencies exceeding those of the 
current. 
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kx (rad m 


Figure 5. Power spectral density of current density Pj{kx) for 
the solar wind interval in Figure 2a of IChen et~ 'all lIMIl') cal¬ 
culated from the trace magnetic spectrum Pb{^x) as Pj(kx) = 
{{V2k^/ sin 9 bv)/Pb (kx) , where sin^Bv is the angle between 
the mean magnetic field and the solar wind (sampling) direction, 
kx = 27r/sc/nsw is the wavenumber in the sampling direction as¬ 
suming the Taylor hypothesis, /sc is the spacecraft-frame frequency 
and Usw is the solar wind speed. The spectrum integrated between 
the black dashed lines (to exclude windowing effects at small 
and instrumental noise at large gives Sj^ = 2.7 X 10“^® 
m“^, which corresponds to a velocity vq = 9.7 km s“^, a factor of 
3.8 smaller than Vg. 

times larger than its rms value. At present it is un¬ 
clear whether the current intermittency in the solar wind 
is strong enough to account for the inferred mismatch 
between the rms electron drift velocity and the ion- 
acoustic speed. It is therefore hard to conclude whether 
the anomalous resistivity can be essential at subproton 
scales. 

5. INTERACTIONS OF ALFVEN WAVES WITH 
ION-ACOUSTIC AND FAST FLUCTUATIONS 

In a low-beta, /3i <C I, and non-isothermal, Ti T^, 
plasma a significant role in a turbulent cascade can be 
played by the ion-acoustic modes and the fast modes, 
where the density fluctuations are not small. Such 
modes are subject to strong Landau damping in the case 
/3e ^ /3i ~ I, which partly explains rather low level of 
density fluctuations observed in the corresponding tur¬ 
bulence. In the non-isothermal, low-beta case, however, 
the density fluctuations associated with these compress¬ 
ible modes become dynamically essential. One of the im¬ 
portant implications is the possibility of parallel energy 
cascade due to decay interactions of the Alfven waves 
with the ion-acoustic and fast waves. 

5.1. Scatter of Alfven waves by ion-acoustic fluctuations 

For an estimate we consider weak turbulence of shear- 
Alfven and ion-acoustic mod es (excellent accounts of this 
theor y can be found in, e.g., lKuznetso^l2001l : iChandranI 
120081) . In the absence of the ion-acoustic modes, the tur¬ 
bulence of shear-Alfven modes is known to transfer the 
energy in the field-perpendicular direction, populating 
mostly wavemodes with k± ^ fc||. The rate of nonlin¬ 
ear energy transfer due to Alfven wave interactions is 
estimated as 7 ^ v^/k\iVA le.g.. lNg &: Bhattacharied 
[mi), where v is the amplitude of the velocity fluctua¬ 
tions in the wave. This expression is applicable also in 



Figure 6. Histogram of the ratio of normalized magnetic field 
fiuctuations to ion d ensity fluctuations for the low beta solar wind 
intervals discussed in IChen et al.l II2014I1 in the range of spacecraft- 
frame frequencies 3 X 10~® Hz to 1 X 10“^ Hz (within the MHD 
inertial range). Absolute values of density fiuctuations are used. 

the case of strong turbulence, if one assumes that the 
critical balance condition v/va ^ k\\/k± is satisfied (e.g., 
iCIoldreich fc SridhaillI995l ). 

The Alfven wave interaction with the ion-acoustic fluc¬ 
tuations is dominated by scattering of the Alfven waves 
by the density inhomogeneities, as the magnetic fluctu¬ 
ations are small in the low-beta ion acoustic modes. In 
the Elsasser variables, = v±b/\/ 47 rp, the MHD eq ua- 
tions take the form (e.g.. iMarsch fc Mangene^lI987fl : 

dtz^ T {VA ■ V)2± -f (2^ • V)2± = -iVp - 
—VA^ ■ V -f \{z~ — z+)V • (&/v^47rp + v/2) + 

+ ^{z-- z+)V ■ VA, ( 35 ) 

dtp+V ■ pv = 0 , ( 36 ) 

where va = B^j and the pressure p includes both 

the kinetic and magnetic parts. The interaction of the 
shear-Alfven waves with the density fluctuations is de¬ 
scribed to the leading order by the last term on the right 
hand side of (I55|) . If the majority of density fluctuations 
have fc[| ^ fe_L, the energy of Alfven waves is transferred 
in the k space mostly along the field-parallel direction. 
Let us assume that an Alfven wave, say propagates 
along the background magnetic field through a series of 
density humps of typical scale I, having intensities Sni. 

Following a standard weak-turbulence consideration, 
we estimate that in each crossing time ~ I /va the ampli¬ 
tude of the wave is distorted by ^ {Sni/n)z~^. Since 
these ampitude changes are random, one needs to accu¬ 
mulate ^ {n/Snifl distortions to change the amplitude 
significantly. During an interaction the frequency of the 
wave changes by Aui ~ ±ujy/2fl/, with the sign depending 
on the mutual directions of propagation of the Alfven and 
ion-acoustic perturbations. The energy diffusion in the 
phase space is thus impeded by a factor {Auj/ujfl ^ fte- 
The inverse time of energy cascade in the phase space 
can thus be astimated as 75 = l3e{vA/l){5rii/nfl. 

Parenthetically, we note that this consideration may 
need to be modihed if the density fluctuations are allowed 
to steepen and develop shocks. In this case, we may as- 
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sume that the wave is propagating through a set of 
density discontinuities of amplitude 6ni separated by dis¬ 
tance 1. In each crossing time ~ 1/va, the wave looses its 
energy; a fraction of energy transferred to the reflected 

wave is equal to 5{z'^Y/{z'^Y ^ {Sni/nY. The fre¬ 
quencies of the reflected waves differ by Aw Ri ±wi/2^ 
from the frequency of the propagating wave. Since the 
signs of individual frequency shifts are random, one needs 
to accumulate {uj/AujY ~ l/(2/3e) subsequent reflections 
before the frequency changes by its magnitude. There¬ 
fore, the time of the frequency shift is St ^ (//nA)/(2/3e)- 
During this time, a fraction of energy [Sni/nY will cas¬ 
cade to smaller scales. The rate of energy transfer is 

therefore 75 = [5{z+Y/5t\/{z+Y ~ l5e{vA/l){Sni/nY , 
which provides the same estimate as before. 

Comparing the rates of the energy transfer, we get: 

is/lA ^ /3e(fc||/fc_L)^(^n/n)2((56/So)"^ (37) 

where we have substituted I ~ l/fc||- Observations 
(|Chen et al.ll2014fl demonstrate that in the regime 1 » 
Pi, the density and magnetic fluctuations are on the 
same order in the MHD inertial interval, Sn/n ^ Sb/Bo, 
see Fig. [51 therefore the Alfvenic and density fluctua¬ 
tions exhibit mostly field-parallel energy cascade inside 
the cone fc||//cj_ > Ijy/P^. Assuming that the veloc¬ 
ity and magnetic fluctuations have the same turbulence 
spectra, we then obtain from the constant rate of energy 
cascade, e ^ -jz"^ ~ const, the energy spectrum 

i?(fc||)(ifc|| oc fc|| (38) 

This result agree s with the analytic consideration of 
([Kuznetsovl 1200111 . if we assume that the propagation 
is mostly parallel. According to the dispersion relation 
(I 25 I) in this case the cascade will have a break point at 
fc ~ fc|| = 1/di. This may be consistent with the obser¬ 
vations in (jChen et al.ll20l4l . 

5.2. Interactions of Alfven waves with fast modes 

An efficient energy transfer in the parallel direction 
can also be provided by the interaction with the fast 
modes. These modes have comparable normalized mag¬ 
netic and density fluctuations, therefore, the interaction 
is described by the last term in the left-hand side and by 
the last two terms in the right-hand side of (l35ll . The 
energy transfer is isotropic in the k-space, with the rate 
estimated as jf ^ {vA/^){Sn/nY, where A ^ 1/fc. This 
rate is comparable to the rate of field-perpendicular en¬ 
ergy transfer due to the shear-Alfven modes. There¬ 
fore, interactions with fast modes provide an efficient 
energy transfer in the field-parallel direction, and can 
explain the spectral break at fc ^ fc|| ^ 1/di. Tbe en¬ 
ergy spectrum, estimated from the constant energy flux 
'ypz'^ = const, is given by 

E{k)dk oc fc-3/2dfc, (39) 

which is also consistent with the analysis of ([Chandranl 
120051 ICho &: L azariaD||200^ , and may agree with the 
observations of fChen et al.l (1201411 . This consideration 
holds except for the propagation directions close to the 
background magnetic field, where the density and the 
field-parallel magnetic fluctuations of the fast mode be¬ 
come small, Sn/n = Sb^fBo = {k±/k)6b/B q. At the 
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Figure 7. Histogram of the ratio of normalized magnetic field 
strength fluctuations to ion den sity fluctuations fo r the low beta 
solar wind intervals discussed in IChen et al.l II2014I 1 in the range of 
spacecraft-frame frequencies 3x 10“® Hz to 1X 10“^ Hz (within the 
MHD inertial range). Note deviations from the pressure balance 
condition. Lower values of the argument are consistent with the 
presence of ion-acoustic modes, higher values - with the presence 
of fast modes. Absolute values of fluctuations are used. 

propagation angles k±/k Y, \fPe, the energy transfer is 
then dominated by the scattering of the shear-Alfven 
waves by the ion-acoustic modes, as described in the pre¬ 
vious section. 

6. OBSERVATIONAL EVIDENCE FOR DYNAMICALLY 
IMPORTANT DENSITY FLUCTUATIONS 

The measurements of spectral anisotropy a re not avail¬ 
able for the low-beta intervals analyzed in (jChen et al.l 
l2014|) . We can however infer the signihcance of den¬ 
sity fluctuations in the energy cascade in a low-beta 
plasma from comparison of the density fluctuations with 
the fluctuations of magnetic field strength. If the den¬ 
sity fluctuations are not dynamically important, e.g., 
they are p assively advected by the shey-Al^ en turbu¬ 
lence (e.g.. lChandraDll20?i^lHowes et al.ll2012t) . they sat¬ 
isfy locally the pressure-balance condition (5|i?|/i?o = 
— (/3e/2)i5n/no- In this case, the fluctuations of the mag¬ 
netic field are dominated by the field-perpendicular com¬ 
ponent, 5bF ^ 5bz, while the small field-parallel fluc¬ 
tuations are related to the fluctuations of the magnetic 
field strength 5\B\ k, 5b z- The observations however 
show that the pressure-balance condition is not gener¬ 
ally observed in the considered case, see Fig. [71 The 
deviations from the pressure balance may be explained 
by the presence of the fast modes, for which 5\B\/Bq = 
Sn/uf), and/or almost field-parallel ion-acoustic modes, 
for which 5\B\/Bej = —{k±/kY{Pe/2)5n/no, which lends 
support to the importance of these modes in the turbu¬ 
lent cascade, as discussed in the previous sections. 

7. CONCLUSIONS 

We propose that according to the critical balance con¬ 
dition, the breaks in the turbulent energy spectra of a 
weakly collisional plasma, as those observed in the so¬ 
lar w i nd and in the inte r stella r plasma (e.g.. [Bale et ahl 
20051: iBourouaine et al.l l2()12t [Markovskii et al.l 120081 
Haverkorn fc Spanglerll2013f) . are related to the breaks in 
the corresponding dispersion relations. For a high beta 











































plasma, /3i ^ 1, the break for the shear-Alfven mode oc¬ 
curs at scales fcj_ ^ 1/pi- For a low beta plasma ,5^ <C 1, 
on the other hand, the situation crucially depends on the 
obliquity of propagation, and the electron plasma beta 
/Se¬ 
lf kl/k^ < /3e, the turbulent energy cascade is predom¬ 
inantly in the field-perpendicular direction. The energy 
transfer is dominated by interactions among the shear- 
Alfven waves, and the spectral break for the shear-Alfven 
mode depends on the combination of the ion-acoustic 
scale and the ion inertial length scale, as follows from 
formula (051) . For /3e < 1 and kz ^ the the energy 
cascade governed by interactions of shear-Alfven modes 
with fast modes can compete with the cascade governed 
by mutual shear-Alfven interactions. In this case sig¬ 
nificant fractions of energy can cascade in the direction 
of large ku. Finally, for k1/k\ > l//3e the energy cas¬ 
cade may be predominantly in the direction of large /cj|, 
and it may be dominated by interactions of shear-Alfven 
waves with the ion-acoustic modes. In the last two cases, 
a spectral break would appear at the scale k ^ l/d,;, 
in ag reement with the recent observations (jChen et al.l 

[Ml . 

Our analysis, together with the observations of 
(|Chen et al.ll20ll . may provide an evidence for the sig¬ 
nificance of a field-parallel turbulent cascade in a low 
beta plasma. They also emphasize the importance of dis¬ 
tinguishing among various plasma regimes in describing 
turbulence in a collisionless plasma, and in interpreting 
observational results. 
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